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The statistical mechanics described by partition lunctions for quantum gravity that assume 
a finite number of degrees of freedom for a fixed finite spacetime region are discussed. It 
is shown that a cosmological constant arises naturally as a Lagrange multiplier needed to 
fix the 4-volume in the canonical ensemble and that in the grand canonical ensemble, the 
cosmological constant naturally emerges in terms of a volume normalized chemical potential 
as A = — 87r/i. This result is demonstrated classically by defining the free energy in terms of 
the classical action and calculating the change in free energy when a unit volume is added 
to a spacetime region having a fixed boundary. 



I. INTRODUCTION 



The cosmological constant is often interpreted to be a result of vacuum energy fTHS] and, as is well 
known, this leads to estimates of the cosmological constant many orders of magnitude larger than 
observations allow. The problem with interpreting the cosmological constant as vacuum energy has 
led to the introduction of alternative models j3HS] such as quintessence jjj. We will argue that the 
notion that the the cosmological constant is determined by vacuum energy rests on an incorrect 
assumption regarding how Einstein's equation arises as a classical approximation to an underlying 
quantum theory. In particular, Einstein's equation follows as a stationary phase approximation to 
a quantum mechanical partition function. If there are an unlimited number of degrees of freedom 
in any region of spacetime, variation of the metric in the action is unrestricted and the resulting 
Einstein equation contains a stress-energy tensor that includes an absolute energy density. The 
zero-point energy would then appear to contribute to the stress-energy tensor. However, if we 
assume that there is only a finite number of degrees of freedom in any finite region of spacetime, 
the variation of the metric must be restricted. We will show that implementing such a restriction 
introduces an energy density scale so that the energy density appearing in the resulting Einstein 
equation can be regarded as a relative energy density that does not include vacuum energy. 

In [8], we derived a partition function for quantum gravity where any finite spacetime region 
contains only a finite number of degrees of freedom. Here we study the consequences of the general 
structure of this theory for the cosmological constant. In Section [TTJ, we show that the cosmological 
constant naturally arises from the canonical ensemble as a Lagrange multiplier, A, needed to enforce 
a fixed 4-volume constraint. That the cosmological constant could be a Lagrange multiplier was 
considered in [6] where it was argued that the approach is problematic because of the need to 
restrict the range of the Lorentzian time integral. Such difficulties do not arise in our approach, 
however, where we work in a local Euclidean spacetime region [Sj. In Section III we discuss the 
grand canonical ensemble and relate the cosmological constant to a chemical potential. We then 
show in Section IV that the conclusion from the previous sections that A = —fi = A/(8tt) follows 
at the classical level by calculating the change in free energy when the 4-volume in a region with 
a fixed boundary is varied. We discuss some of the implications of our results in Section |V| where 
we argue that zero-point energy plays no role in determining the cosmological constant. 



* This paper is available at 



http://archive.org/details/CosmologicalConstantChemicalPotential 
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II. PARTITION FUNCTION FOR QUANTUM GRAVITY 

In [H] , we derived an expression for the canonical ensemble partition function that can be written 
in the form 

Q(V,E) = tr e~ s (1) 

where V is the fixed 4-volume in a region M bounded by the 3-surface £ = dM and the trace 
involves integration over a finite number of metric, ghost, and matter degrees of freedom. The 
4-volume was required to be held fixed, because it was assumed that there is a finite number of 
degrees of freedom for any given finite 4-volume and fixing the 4-volume allowed the partition 
function to be expressed in terms of a finite number of integrations. Classical general relativity 
results from taking the stationary phase approximation to the partition function with a Lagrange 
multiplier A introduced to keep the volume fixed. The Euclidean spacetime action can be written 

S = S G + S M (2) 

where Sq is the gravitational action, which may be taken to a a suitably Wick-rotated version 
of the Einstein-Hilbert action [SI |5] and Sm is the action for matter fields. For considering the 
stationary phase approximation, we may take (in d = 4) 

S G = -^rl d 4 x^R-±- [ d 3 xVhK = S B + S s (3) 
16vr J M v 8vr J 9M 

where R is the Euclidean Ricci scalar, h is the determinant of the metric induced on dM, K is the 
extrinsic curvature of dM, and Sb and Ss refer to bulk and surface terms respectively. We may 
write the matter contribution as 

S M = - j d 4 x^/g C m (4) 

Defining 



(5) 



2 5(^gC m ) 
\fg 5g ab 

and imposing the stationary phase condition 

J-{S G + S M + XV) =0 (6) 
og a b 

results in the Einstein's equation 

G ab = 8iT(T ab - Xg ab ) (7) 

so that 87rA is the cosmological constant. We will show in the next section that the cosmological 
constant is related to a suitably normalized chemical potential. 

III. GRAND CANONICAL PARTITION FUNCTION 

Here we consider a grand canonical ensemble partition function. Before describing the thermo- 
dynamics, we first discuss the geometry of the statistical mechanical system. As we discussed in 
[8J, a Lorentzian metric tensor can always locally be written as 

d 2 s = -dt 2 (l + n(x)) + (5ij + g ij (x))dx i dx j + v^dx^t. (8) 
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where n(x), gij(x) and V{(x) are quadratic in spatial coordinates and are independent of the time 
coordinate t and cubic and higher terms have been dropped. This can be understood as the Wick 
rotation of the Euclidean metric 

d 2 s = dr 2 (l + n(x)) + (<% + g ij (x))dx t dx j + Vi{x)dx l dT. (9) 

where Wick rotation takes t — > —it and Vi — > ivi. The r coordinate is naturally interpreted as the 
Euclidean time parameter in a thermodynamic description of spacetime degrees of freedom with 
r running from to the inverse temperature (3. It is convenient to define S as follows. Define a 
2-surface So a t a time slice r = and extend this surface along r to form a cylinder. S may be 
defined as the 3-surface bounding the 4-volume contained within the extended surface So and the 
surfaces r = and r = (3. 

To define the grand canonical ensemble partition function, we start with the canonical ensemble 
partition function and sum over discrete possible 4- volumes, V, with the partition function weighted 
by a factor of e^ v where n is the chemical potential with a volume normalization. The grand 
canonical partition function can then be written as 

Z(fi, 0, S ) = tr e~ SE+llV (10) 

where the trace now include a sum over 4-volume and the number of terms included in the expansion 
of the metric for a given 4-volume depends on the 4-volume. The stationary phase approximation 
then gives Einstein's equation with A = — 87r/i. 

We will now show that the chemical potential [i appearing in this partition function can be 
related to a partial derivative of the free energy by taking the differential of In Z(/x, /?, So) and 
comparing to the first law of thermodynamics. Even though this analysis is very similar to the 
traditional case|10|. we will work through this is some detail since our normalizations differ from 
that of standard statistical mechanics. 

We can define a thermodynamic average of the energy by 

£ = -^lnZ(/x,/3,S ) (11) 



and similarly it follows from equation (10) that the average 4-volume is 



V = ^lnZ(/x,/3,S ). (12) 

The pressure in conventional statistical mechanics would be defined as P = |j d y( 3 ) InZ. Here, the 
boundary Sq plays the role of volume in ordinary statistical mechanics. We can define a pressure 



as 



/3 3F So 

where the volume derivative q^t is defined by 



P =\l^^ Z ^P^) (13) 



-0 



e f&tjm (14) 

where Si is the surface So with each point moved a uniform infinitesimal distance normal to So 
such that the 3-volume between Si and So is dV = Vsi — Ve q . The differential of In Z(fi, (3, So) is 
then 

din Z = -Ed/3 + (3PdV So + Vdfi (15) 
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from which we obtain 

dE = -d(ln Z + PE-pV)- PdV^ + -fidV. (16) 

Writing T = 1//3 and S = In Z + f3E - p,V, we have 

dE = TdS - PdVs + TpdV. (17) 
This can be recognized as the first law of thermodynamics with a modification to the normalization 



of the chemical potential term. It follows from equation (17) that the chemical potential can be 



written in terms of the Helmholtz free energy, A = E — TS, as 

where we have dropped the bar on V for simplicity of notation and it should be understood that 
the thermodynamically averaged V is used here. We see that p is fH times the change in free energy 
when a unit of 4- volume is added to the volume bounded by the surface S holding the boundary 
and temperature fixed. If we define a 3-volume as = V//3 we can write 

and we see that p can be described as the change in free energy when a unit of 3-volume is added 
to the volume bounded by the 2-dimensional surface So holding the boundary and temperature 
fixed. This is naturally interpreted as a volume-normalized chemical potential. Expressing the free 
energy in terms of the canonical ensemble partition function, we find 

M = -^lnQ(F,E). (20) 

Since we have found A = —8irp and A = 87rA, it follows that 

A = ^lnQ(V,E) = -/i. (21) 

In the next section, we will calculate the partial derivative of the canonical ensemble partition 



function at the classical level and show that equation (21) holds 



IV. CHEMICAL POTENTIAL IN THE CANONICAL ENSEMBLE 

Consider a slowly varying, spherically symmetric distribution of matter that can be described 
as an ideal fluid with density p and pressure p. We will calculate the free energy in a finite volume 
region where the length scale of the volume is large compared to but small compared to the 
length scale on which p and p vary. We will hold the metric, p and p fixed on the boundary of 
our volume and vary the volume inside the boundary by varying the metric. We do this at the 
stationary phase level of approximation for the free energy and calculate the chemical potential 



using equation (20) 



In keeping with the spirit of the local thermodynamic equilibrium approximation, we assume 
that the length scale, Ly associated with our fixed volume region is small compared to the curvature 
length scale L c . If we were to calculate ^ lnQ(V, S) as a power series in p, p and A only the linear 
terms would contribute in the limit that Ly « L c since higher order terms contain factors of 
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Ly /L c . The derivative with respect to V, however, can remove a factor of A from ^ lnQ(V, S) so 
to guarantee that we get all linear terms in ^\nQ{V, £), we will keep terms to quadratic order 
in In Q(V,£). 

We cannot assume that p and p are constant since we would not be able to satisfy Einsteins 
equation if this were the case. In particular, the Tolman-Oppenheimer-Volkoff equation|llj. which 
follows from the Einstein equation, requires that the pressure depend quadratically on radial coor- 
dinate r for small f for a slowly varying spherically symmetric distribution of matter. If we assume 
that the matter obeys an equation of state such that the pressure is a function of energy density 
then the energy density must also vary quadratically with f. We may therefore write 

P = Po ~ Pif 2 = Pb ~ Pi{r 2 - ff) (22) 

where f b is the value of f at the boundary, po , p b are the densities at r = and at the boundary 
respectively, and pi is a parameter describing the variation of p. Similarly, we may write 

P = Po -Pir 2 = Pb -pi(r 2 - rf). (23) 
From the Tolman-Oppenheimer-Volkoff equation, we have to quadratic order 

1 2 

pi = 2vr(p b + p b ) (p b + -p b - -A) (24) 
The density variation parameter can be related to p\ by the equation of state 

« = i£k (25) 

where dp/ dp can be taken to be a constant in evaluating p\ to quadratic order. 
The solution to Einstein's equation can be written as[llj 

ds 2 = e 2<t> dT 2 + ( (1 _ ~ l d ~2 + -2 d0 2 (26) 

where 



dcj) m(f) + 4-Kr 3 (p — A) 
dr f(r — 2m(f)) 



(27) 



m(r) 



4vr / (p(r) + \)r' 2 dr' (28) 
Jo 



and r = (A/Att) 1 ^ 2 where A is the area of the 2-sphere of radius f 

We consider an unperturbed geometry having some 4- volume V(X) and consider a perturbed 
geometry with a 4- volume V(X + 5X). For both geometries, the metrics can be written in the 



form of equation (26). A radial coordinate f = (A/ Ait) 1 / 2 can be defined for each geometry, 
however we cannot simply identify coordinates in the two geometries except at the boundary where 
the geometry is held fixed. For the angular components, matching the metric on the boundary 
is achieved simply by the condition that f = r b at the boundary in both cases. For the time 



component, this is achieved by imposing a common boundary condition in equation (27). For the 
radial component, we cannot set <?f,f in the unperturbed geometry equal to g^f in the perturbed 
geometry on the boundary. This is because gff measures distances off of the boundary and we 
cannot, in general, identify points off of the boundary labeled by the coordinate f in the two 
different geometries. However, the choice of coordinates is arbitrary and we need not express 
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the metric in both geometries using the f = (A/ An) 1 / 2 coordinate. Let write the metric in the 
unperturbed geometry using the radial coordinate as r = r and in the unperturbed geometry we 
write the metric using a coordinate r 1 defined by 



f = f(r') = r + 7(7*' — r b )(r' /r^) 



(29) 



Note that on the boundary f = (A/ An) 1 / 2 = r' which can be identified with r = (A/ Air) 1 / 2 on the 
boundary since the metric and hence the area A is held fixed on the boundary. We can choose 7 
so that the metric expressed in terms of r' is equal to the metric expressed in terms of r on the 
boundary, i.e., g rr = g r > r i on the boundary. We can then identify r' with r up to an arbitrary 
diffeomorphism that preserves the boundary and the metric on the boundary. Since the theory is 
invariant under such boundary preserving diffeomorphisms, we will henceforth drop the distinction 
between r' and r. 



Using equations (22) and (23) in equations (26) - (28), the metric to second order can be written 



as 



ds 2 = c(l + af 2 + -(a 2 + b)r 4 )dr 2 + {dr/dr) 2 {l + Ar 2 + Bf 4 )dr 2 + f 2 dn 2 (30) 



where f = r for the unperturbed geometry, r is given by equation ( 29 ) with r' replaced with r for 
the perturbed geometry, 



Att 

a = —(Pb + A + pir b ) + A7T(p b - A +p±r 2 ) 



(31) 



4"7T 47T 9 

— pi - Air Pl + 2(— ) ( Pb + A)(p 6 + 3p 6 - 2A), 



(32) 



.4 



y(p 6 + A + pir b ), 



(33) 



_ 87T / 8tT ■. 2 , . s ' 

b = -— Pi + (—) (p b + xy 



(34) 



and c is a constant that is determined by the boundary conditions. We will denote values for the 
unperturbed geometry with a subscript. For the unperturbed geometry, we take cq = 1. The 
value of c in the perturbed geometry is then given by requiring the metric to remain fixed on the 
boundary. c 1//2 , which will be need shortly, can be expresses to second order as 



.1/2 



1 + 2^0 - a)r 2 b + ^(a - a) 2 r fe 4 + ^(b - b)4 



Finally, 7 is determined by matching g rr at the boundary. We find 

'2 1 73 J2 



1 



1 (( l + Apij + Bo 



2r b \\ l + Ar 2 + Br 2 



1 



(35) 



(36) 



Since in the stationary phase approximation (in our sign convention) Q = e 



(37) 
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and we can calculate the chemical potential by taking the volume derivative of the total action 
evaluated at the stationary phase condition. We calculate three contributions to the chemical 
potential. The first is the contribution from the volume term in the action. This can be calculated 
using the trace of Einsteins equation 



i? = 32vrA + 8vr( / 9-3p). 



(38) 



The contribution to the free energy from this term can be obtained by substituting the previous 
expressions for p and p and using the volume element determined from the above metric. The 
derivatives with respect to volume can be written as 







1 







(39) 



dV dV/dX dX 

After some algebra, we find that the contribution to the chemical potential from the bulk term is 

3 



Sb 



33 107 9 (p b +p b ) 
+ 2X -V0 Pb + ^0 Pb -35^p7dp- 



(40) 



There is also a contribution to the free energy from the term from the surface part of the action. 
This can be calculated from 



d 3 xVhK 
dM g rr 



1 d 



1/2 dr -r\r=r b 



(41) 



where A r is determined from the metric. Using equation (30) in equation (41) and applying 
equation(39), we obtain (again, after some algebra) 

3 



Ss 



2vrrf 



3A + -p 6 



27 9 jpb+Pb) 
U Pb 35 dp I dp ' 



(42) 



Finally we have a contribution from the matter part of the action. This can be derived from 
equation §5§. We have 



d_ 
dX 



S 



M 



d_ 
dX 



(43) 



Using equation (39) and using equation(30) in equation(43) with the stress-energy tensor of an 
ideal fluid, we find that the contribution of the matter part of the action to the chemical potential 
of the spacetime degrees of freedom can be written as 



Sm 



Adding these three terms together we have simply 



p = p S B + ^Ss + ^Sm 



A 
'8tt' 



(44) 



(45) 



V. COSMOLOGICAL CONSTANT 

A common interpretation of the cosmological constant is that it describes the energy density 
of the vacuum. From the point of view of the present work, this is a misconception that arises 
from the incorrect notion that general relativity depends on an absolute energy density rather 



s 



than a relative energy density. As noted earlier, this is consistent if there are an unlimited number 
of degrees of freedom in any finite region of spacetime. In our approach, however, this is not 
the case and a relative energy enters into general relativity just like every other area of physics. 
In particular, it is the energy density relative to a chemical potential that is relevant in general 
relativity. 

To explore this further, consider an analogy with condensed matter physics. For non-interacting 
bosons, the occupation number can be written as 

nk = e/^-ri-i ' (46) 

where are the energy levels of the system. If we don't include zero-point energy in the definition 
of €k then we should not include it in the definition of /i. The zero-point energy represents an 
arbitrary zero of energy and as long as and \x are defined relative to the same choice for the 
zero of energy, the occupation number does not depend on this choice. The same is true in general 
relativity. In synchronous coordinates in Lorentzian spacetime (<?oo = — l)j we have 

Goo = 8^r 00 + A = 8^(T 00 - fi) . (47) 

We see that as long as Too does not include zero-point energy then neither should \x. If one insists 
on adding a constant energy density to Too to represent zero-point energy, one should add the 
same constant to \x so that Too an d [i are defined relative to the same arbitrary zero of energy. But 
then this constant energy density does not enter into Einstein's equation. Since the conventional 
definition of Too m general relativity does not include zero-point energy, then neither should /U. More 
generally, it is only the difference between T a ^ and —g a b H which enters Einstein's equation and we 
should either include zero-point contributions to both T a b and —g a b H or to neither. Problems only 
occur if we attempt to include zero-point contributions in one of these terms and not the other. 

A general feature of discrete spacetime theories [8] is the existence of a finite number of degrees 
of freedom associated with any given finite region of spacetime and as we have seen a cosmological 
constant that is unrelated to zero point energy enters into such theories. The observed evidence for 
a cosmological constant with a magnitude many orders of magnitude below the Planck scale can be 
taken as an indication that discrete degrees of freedom underlie spacetime. Note that although our 
approach divorces zero-point energy from the cosmological constant and thus explains why naive 
estimates are off by many orders of magnitude, it does not explain why the cosmological constant 
takes on the particular value that it does and it does not address the coincidence problem, though 
it is not clear how serious a problem this is|12j. 



VI. CONCLUSION 

We have shown that if one assumes that there should be a finite number of degrees of freedom 
in any finite spacetime region, then the existence of a cosmological constant follows quite generally 
from the partition function written as a sum over the degrees of freedom in a fixed 4- volume where 
a Lagrange multiplier must be introduced to fix the 4-volume when the metric is varied. The 
cosmological constant written in terms of the Lagrange multiplier is simply A = 87rA. 

If we assume that a thermodynamic interpretation of the finite number of degrees of freedom is 
appropriate, as we argued in jS], then the cosmological constant is naturally interpreted in terms of 
a volume-normalized chemical potential. We showed that A = — 87r/x. This relation follows directly 
from the grand canonical partition function. We explicitly verified that fi = —A for a classical, 
slowly- varying, spherically symmetric distribution of matter by calculating InQ in the stationary 
phase approximation and differentiating with respect to 4-volume to obtain the chemical potential. 
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We discussed the implication of the interpretation of the cosmological constant as a chemical 
potential. In particular, we argued that this implies that there is no zero-point contribution to the 
cosmological constant, contrary to a popular misconception. The theory presented here does not 
shed any light on the particular value that the cosmological constant takes other than to show that 
estimates that are off by 60 or 120 orders of magnitude are based on an incorrect understanding 
of how zero-point energy enters, or fails to enter, Einstein's equation. From the point of view 
presented here, the cosmological constant can be interpreted as a parameter that describes how 
much 4-volume in empty spacetime is contained inside a region bounded by some fixed boundary 
(i.e., how flat empty spacetime is), just as the chemical potential for particles in ordinary statistical 
mechanics is a parameter that describes how many particles are within a fixed region. 
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